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Abstract. In classical Drude theory the conductivity is determined by the mass of the 
propagating particles and the mean free path between two scattering events. For a quantum 
5_l ■ particle this simple picture of diffusive transport loses relevance if strong correlations dominate 

' the particle motion. We study a situation where the propagation of a fermionic particle is 

possible only through creation and annihilation of local bosonic excitations. This correlated 
■ quantum transport process is outside the Drude picture, since one cannot distinguish between 

free propagation and intermittent scattering. The characterization of transport is possible using 
the Drude weight obtained from the f-sum rule, although its interpretation in terms of free mass 
' O ' and mean free path breaks down. For the situation studied we calculate the Green's function 

and Drude weight using a Green's functions expansion technique, and discuss their physical 
O . meaning. 

(N 
> . 

^ . 1. Motivation 

\ The motion of a fermionic particle, such as an electron or hole, that interacts strongly with some 

■^J- ■ background medium is a constantly recurring theme in condensed matter physics. The present 

\q ■ contribution is motivated partly by the study of the motion of a hole in an antiferromagnetic 

spin background, as in a doped Mott insulator. Mott insulators occur for electrons on a lattice 
subject to strong Coulomb repulsion p], and are commonly studied in the context of the Hubbard 
model. Since putting two electrons at the same lattice site costs a large energy if Coulomb 
repulsion is strong, the number of doubly occupied lattice sites is negligibly small. At half-filling 
^ . we can therefore assume that every lattice site is occupied by a single electron. The system 

is an insulator because electron motion, which creates doubly occupied sites, is energetically 
forbidden. 

If electron motion is suppressed, the electron spin is the only relevant degree of freedom. 
Perturbation theory in the kinetic energy shows that spins couple antiferromagnetically in a Mott 
insulator. Starting from the Hubbard model, spin dynamics is described by the antiferromagnetic 
Heisenberg model with complex properties arising from the interplay of strong correlations and 
quantum spin fluctuations. In a very simple picture, we may alternatively consider the Neel 
state of a classical antiferromagnet (see Fig. Q] upper panel). 

If an electron is removed from a two-dimensional (2D) Mott insulator in a Neel state, 
the hole can move without creating doubly occupied sites. But the moving hole distorts the 
antiferromagnetic spin background, creating misaligned spins with substantial energy (see Fig. [TJ 
lower panel). The 'string' of misaligned spins that forms along the path of the hole strongly 
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Figure 1. Upper panel: Hole in a Neel antiferromagnet on a 
square lattice. Lower panel: When the hole moves by nearest- 
neighbour hopping, it creates misaligned spins that do not 
match the Neel order. These distortions of the spin background 
lead to parallel spins on neighbouring sites - indicated with a 
dashed line - with a substantial increase in energy. To 'unwind' 
the string of background distortions, the hole in principle has 
to retrace the path it moved, and is therefore bound to its 
origin. The existence of (quantum) spin fluctuations, and spin- 
order restoring processes (see Fig. |2J) weakens, or completely 
destroys, this 'string effect'. 



restricts propagation, and tends to bind the hole to its origin. Therefore the hole does not move 
as a free particle. In a quantum antiferromagnet, which does not possess the strict order of the 
Neel antiferromagnet, the spin distortions can 'heal out' or relax by quantum spin fluctuations. 
It is then the relaxation rate of spin distortions that determines how fast the hole can move: The 
hole is continuously creating a string of distortions but can move on 'freely' at a speed which 
gives the distortions time to decay. 

In the absence of quantum spin fluctuations free-particle like motion of a hole is completely 
suppressed. Surprisingly enough, a hole can propagate through so-called 'Trugman paths' [2], 
which realize a certain combinatorial rearrangement of spins compatible with the spin order. The 
simplest of these processes, where the hole moves by two sites and restores the antiferromagnetic 
spin alignment, is shown in Fig. [2j In contrast to the situation for an almost free particle 
propagation through such processes depends entirely on the successive rearrangement of spins 
in precise coordination with the hole motion, just like executing a complicated dance pattern. 
With respect to the correlations between spin and hole motion we call this a correlated transport 
process. 

The focus of the present contribution is on similar correlated transport processes in a simpler 
fermion-boson model. This model is introduced in the next section (Sec. [2}. In Sec. Owe define 
the Drude weight as the appropriate measure of the efficiency of a correlated transport process. 
In Sec. H] we derive the one-fermion Green's function, which is discussed together with the Drude 
weight in Sec. [5j before we conclude in Sec. [H 

2. The model 

The picture sketched here — how the motion of a particle depends on the creation of background 
distortions and their decay — is very general, and applies e.g. to the polaron problem of electron 
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Figure 2. The simplest 'Trugman path' that allows for propagation of a hole in a Neel 
antiferromagnet [2]. After six steps, the hole has moved by two sites to a next-nearest neighbour 
and the Neel spin order is restored. 
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Figure 3. Lowest order vacuum-restoring process that allows for propagation of a fermion in the 
model Eq. ona ID chain (upper row) and 2D square lattice (lower row). In any dimension, 
it comprises six steps and propagates the fermion by two sites. In steps 1-3, three bosons are 
excited, which are consumed in steps 4-6. The efficiency of this process consequently scales as 
tl/ul (see Sec. I5TT]) . Note that the 2D process is the exact translation of the Trugman path 
shown in Fig. [21 where the correlated many-particle background of the spin model is replaced 
by the bosonic vacuum. For our model, such processes play a role already in ID. 

motion in a deformable lattice studied using the Holstein model [3,4J, and to charge transport 
in the presence of strong electronic correlations, such as for hole-doped Mott insulators studied 
using the t-J-model [3(6]. For that reason, we recently introduced a simple but general fermion- 
boson model that captures the fundamental aspects without modelling any very specific details. 
This model was originally proposed by Edwards [7J and subsequently used to study single-particle 
quantum transport [8] and metal- insulator transitions [9tillj. 
We here consider the Hamiltonian 

* = c j M b l + h) + ^o£ h\b t (i) 

for a spinless fermion (created by cj) coupled to bosons (created by b\) of energy ojq. The first 
term in the Hamiltonian specifies how bosons are created (and annihilated) when the fermion 
hops to a nearest neighbour. Any hop of the fermion changes the number of bosons by one. In 
this sense, fermion motion is boson-controlled. 

Fermion hopping creates a boson only on the site the fermion leaves, and destroys a boson 
only on the site the fermion enters. As a consequence motion of the fermion creates, in a 
lattice of arbitrary dimension, a 'string' of bosonic excitations: In each hop the fermion creates 
additional bosons that can be destroyed only by returning to sites already visited. This 'string 
effect' restricts propagation just as for the hole in the 2D antiferromagnet. 

The original model (in Refs. [ZHS]) contains an additional term — + &i) that accounts 

for boson relaxation. Bosons can decay when A ^ 0, so that in this case the string effect is 
weakened. In the present situation, the string effect is not overcome by boson relaxation. But 
just as for the Trugman paths, there exist processes whereby the fermion propagates, while the 
boson vacuum is restored [5J. The simplest of these vacuum-restoring processes comprises six 
steps, promoting the fermion by two sites (see Fig. [3]). If we compare to Fig. [2j we see that the 2D 
process in our model is an exact representation of a Trugman path. Here, such processes occur 
already in ID, and can be considered as a projection of the corresponding 2D process. In our 
model such correlated transport processes emerge quite naturally from some basic assumptions. 
Our general model now allows us to study the implications of the existence of such processes for 
the motion of a fermionic particle. An analytical calculation of the one-fermion Green's function, 
which is only possible in the A = 0-case considered here, enables us to characterize correlated 
transport using the Drude weight. We restrict ourselves to the study of a single fermion in ID 
at zero temperature. 
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\ gifc" c t | vac ^ = e ~ ik |° Figure 4. Examples of the graphical notation 

^— u n n+ , fermion-boson states. An empty site is 

represented by a dot (•), the fermion by a circle 

} j e lkn c\ l+1 b ] n \\a,c) = e~ lk \* °) (o), and a boson by a star (★). The origin is 

n indicated with an underscore (_). All states 

^ ^ ^ ^ ^ are Bloch states with momentum k. The last 

/ e c n^n-i^n^n+i I vac ) row shows the symmetric state which appears 

n in the 6-step-process (cf. Fig. [5]). 



3. Drude weight 

Transport of almost free particles is conveniently characterized through the particle mass and 
an effective mean free path. How, then, can transport be characterized in a situation where 
it differs entirely from free particle motion? As stressed by Kohn [12] one should consider the 
Drude weight D, which can be defined through its relation 

"-VP".'), (2) 
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to the quasiparticle mass m*, where E(k) is the quasiparticle dispersion. Alternatively, the 
Drude weight can be introduced through, the relation Re — ^ttD S(oj^ H- Creg 

(w) for the real 

part of the conductivity ct(uj), or be obtained from the f-sum rule, which relates the integrated 
conductivity to the kinetic energy. Note that for our model with boson-controlled hopping, the 
current operator is more complicated than usual [8]. 

The Drude weight generalizes the original Drude theory with its concept of an effective 
mean free path between two scattering events, which is valid only for an almost free particle. 
Note that just through the above relation between Drude weight and quasiparticle mass the 
latter is identified as a relevant quantity in a rigorous quantum transport theory. Indeed, the 
quasiparticle mass itself is a rather complicated quantity, which can be strongly renormalized 
as a consequence of correlated many-particle dynamics. 



4. Derivation of the Green's function with a 3-boson constraint 

Even for the simple-looking Hamiltonian in Eq. ([T|) an exact calculation of Green's functions 
is not possible. The situation changes if an additional constraint is imposed, namely a bound 
on the number of simultaneously excited bosons. Since the vacuum restoring 6-step-process in 
Fig. [3] involves 3 bosons, a non-trivial result is obtained allowing for at most 3 bosons in the 
system (we will justify this restriction further in the next section). Then, only a finite number 
of states contribute to the one-fermion Green's function 

G(k,z) = (vac\ck[z — i/] -1 cjjvac) . (3) 



This enables us to derive a closed expression for this Green's function. 




Figure 5. Graph structure for the calculation of G(k,z). Each unlabelled edge corresponds to 
a matrix element tb, the labelled edges to t& times the given factor. The v2 arises from matrix 
elements of bosonic operators. 



4-1- The tree of states 

As in Fig. [31 we use a graphical notation for states, explained in Fig. HI The calculation requires 
the use of Bloch states with given momentum k. As an example, |o) = ^ n e lfcn Cn|vac). The 
graph in Fig. [5] shows all states with at most 3 bosons that are accessible from |o) by repeated 
application of H, while the edges of the graph correspond to non-zero matrix elements of H. 

The graph of states contains a loop corresponding to the 6-step-process. To obtain a 
simple expression for the Green's function with a recursive calculation we must get rid of this 
loop. To this end we perform a change of basis that exploits the reflection symmetry of the 
Hamiltonian ([1]), and rearrange all states in the tree structure shown in Fig. [6) Note that 
the states in this tree are not eigenstates of the parity (or reflection) operator, since they all 
have finite momentum k. Note further that it is important to choose the position of the origin 
correctly. The origin in the two states below the tree root is shifted by one site. Then the 
symmetric state in the last line of Fig. [H which appears at the bottom of the loop in Fig. 
appears only in the left subtree. The corresponding linear combination with a minus sign, 
which would appear in the right subtree, is zero due to reflection symmetry. Therefore Fig. [6] 
has no loop. The shift of the origin accounts for the two sites the fermion propagates in the 
6-step-process. 

4-2. Recursive calculation of the Green's function 

The tree in Fig. fallows for the recursive calculation of G(k,z). The basic elements in this 
calculation are Green's functions G s {z) associated to a state \s) and the subtree with \s) at its 
top. Such a Green's function is defined through the usual equation G s (z) = (s\[z — H s }~ 1 \s) 
(similar to Eq. (J3])), where H s is the projection of the Hamiltonian onto the subspace spanned by 
all states in the subtree. Note that these projected Hamiltonians commute for disjoint subtrees, 
so that the corresponding Green's functions can be calculated independently. It is for this reason 
that we had to remove the loop in the original graph of states (Fig. [5]). 

At each step of the calculation, a Green's function G s (z) associated to a subtree is written 
in the form 

G s (z) = [z-a- bfGdz) - b 2 r G r (z)}- 1 , (4) 

where Gi(z) and G r {z) are the Green's functions associated to the left and right subtrees one 
step further down the tree, bi and b r are the matrix elements of H along the corresponding edges, 



Figure 6. Rearranged tree structure for the calculation of G(k, z) after a basis change. Note 
that the matrix element \/2i& at the edge in the D-subtree results from the scalar product of 
the two corresponding states, and not from bosonic operators. 



and a = (s\H\s) is the diagonal matrix element of H in the state \s). Note that when \s) is at 
the bottom of the tree, the associated Green's function is simply given by G s (z) = [z — a] -1 . 
To justify Eq. (J3|) we proceed in a standard way and split H s as 



H s = a\s)(s\ + b t Vi + Hi + b r V r + H r 



(5) 



where \l), \r) correspond to the states associated with Gi(z), G r (z), and Vi = \s)(l\ + \l){s\, 
V r = \s)(r\ + Multiplying both sides of the operator identity 



l = [z-H r - Hi]- l [z -H s + a\s)(s\ + hVi + b r V r ] 
with [z — Hg]^ 1 from the right, we get the (generalized) Dyson equation 

[z - Hs}' 1 = [z- H r - H^- 1 + [z- H r - H^- 1 [a\s)(s\ + 6,^ + b r V r } [z - Hs}- 1 
Taking scalar products with the state \s) leads to 

G s (z) = z- 1 + z-'faGsiz) + bi(l\[z - Hs}' 1 ^} + b r (r\[z - Hs}- 1 ^} 

and in a second step, with |r), to 

{l\[z-H s }- 1 \s) = Gib l Gs{z) , 
{r\[z-H s }- 1 \s) = G r b T Gs{z) . 



(6) 



(7) 



(8) 



(9) 



Inserting Eq. ([9]) into Eq. ([8]) and solving for G{s) yields our result Eq. (£[]). 

We now apply Eq. (0]) to the tree of states in Fig. [6) Working from the top of the tree of 
states to the bottom, one obtains 



G(k, z)= z- 2t\{K{z) cos 2 k + L{z) sin 2 k) 



z - t 2 b (K(z) + L(z)) - tl(K(z) - L(z)) cos 2k 



(10) 



where the intermediate results to the subtrees labelled in Fig. [6] are 



C(z 
D(z 
E(z 
K(z 
Liz 



z — 2u)q 
z — 2cjo 
z — 2u>o 



z — 3u>o 

3*g 

z — 3u>o 

t 2 
h 

z — 3cjo 



(11) 



z-uj -t 2 b (C{z) + D{z)) 
z-u Jo -t 2 b (C(z)+E(z)) 



4-3. Discussion 

Our derivation of the Green's function is closely related to well-known methods such as the 
renormalized perturbation expansion |13j . the (Lanczos) recursion method [14] or the Mori- 
Zwanzig formalism [15]. For our problem the most important aspect is the appearance of two 
Green's function in Eq. Q, corresponding to the two branches of each node in the binary tree 
of states in Fig. [6J This generalizes the standard recursive construction leading, e.g. in Lanczos 
recursion, to a continued fraction of the form 



G(z) 



1 



a 



(12) 



oi 



a-2 



In the present case this construction would result in very complicated expressions, since arbitrary 
linear combinations of the states shown in Figs. [5l [6] occur in the recursion, e.g. of states with 
different boson number. The additional freedom provided by Eq. is important to obtain 
simple results for G(k, z) in closed form (obviously, our result is not given in the form of 
Eq. ([12]) ). The success depends on our ability to find a suitable arrangement of states in a 
tree as in Fig. EJ This task can hardly be automatized, and such an arrangement is difficult to 
find for more complicated situations. Nevertheless, the idea presented here can be generalized 
so that it may be useful also in further calculations. 



5. Evaluation of the Green's function and Drude weight 

In Fig. [7J we show the spectral function 

A{k,u) = -(l/n)lmG(k,z) . (13) 

We see that the Green's function is periodic in k with period ir. The reason is that the vacuum- 
restoring 6-step process (Fig. [3]) propagates the fermion by two sites, resulting in the term 
ex cos2/c in Eq. ()10p . The spectral function is dominated by a coherent quasiparticle band at 
low energy, which is almost dispersionless and separated by a large gap (of the order ujq) from 
other (coherent) signatures at higher energies. Note that the dispersion of the quasiparticle 
band is related to the Drude weight, as discussed in the next subsection. As Fig. [7] further 
shows, G(k,z) compares quite well to the numerical result without the 3-boson-constraint [8]. 
The reason is that, for all but very small wo/ib — 1, hole motion is dominated by the 6-step 
process that involves only 3 bosons. We will next discuss how this fact shows up in the Drude 
weight. 




CO 



Figure 7. Spectral function A(k,u) for tf, = 1.0 and ojq = 1.0 with the 3-boson constraint from 
Eq. (|10p (left panel) compared to the numerical result (right panel, from Ref. [8]) without the 
constraint. 



5.1. The Drude weight 

Following Eq. ([2]) the Drude weight D is obtained from the ground-state energy 

E(k) = E + (D/2)k 2 + 0{k A ) , (14) 

which is given by lowest pole of G(k, z). With decreasing ujq, the Drude weight grows as boson- 
controlled transport becomes more efficient, but remains finite for ojq — > (see Fig.[8j). Since the 
6-step-process involves 5 intermediate states with excited bosons, the leading order of D/tb in an 
expansion in tb/too is (tb/coo) 5 - Any process contributing to D/tb that has an intermediate state 
with more than 3 bosons, must be of order (ij/cJo) or higher since at least 4 bosons have to be 
excited and de-excited afterwards, which requires 8 steps. Therefore the asymptotic behaviour 
of D for large uio/tb can be obtained from G(k,z) in Eq. (fTUj) . with the result D ~ ^/(3u;q) 
shown in Fig. [SJ 

Once again, we compare our analytical D to the numerical result without the 3-boson 
constraint. Both results fit surprisingly well although the restriction to a small number of bosons 
might appear to be a crude approximation for small wo/tfe- But the inset in Fig. [8J shows that in 
numerical calculations without the 3-boson constraint the number of bosons in the ground-state 
remains small over a large range of values of uio/tb- This indicates, together with the agreement 
between the two curves for D, that the main contribution to fermion propagation comes from 
the 6-step-process even when uo/tb is small, although higher order vacuum-restoring processes 
involving more than three bosons might become energetically favourable in this case. One 
explanation for that behaviour is that correlated transport through vacuum-restoring processes 
relies on the correlated excitation and de-excitation of bosons. Higher order processes involving 
more than 3 bosons cannot give a large contribution to D since correlations between many bosons 
are inevitably reduced. This also explains the saturation of D as ujq — > 0: The contribution of 
the 6-step-process to D remains finite for ooq — > (restricting to 3 bosons, the Hamiltonian is a 
bounded operator, and E{k) and D remain finite). Then, if the main contribution to D comes 
from the 6-step-process, D itself remains finite for loq — > even without the 3-boson constraint. 



Figure 8. Drude weight D with 3-boson 
constraint, and from exact numerics without 
the constraint [8]. The dashed line gives 
the asymptotic result D = tf/(3ujQ) for large 
ujo/tf,. The inset displays the average number 
of bosons in the ground-state without the 3- 
boson constraint. 



6. Conclusions and Outlook 

In this contribution we discuss correlated transport through vacuum-restoring processes in a 
general fermion-boson model. We obtain a simple analytical expression for the one-fermion 
Green's function. ^From the Green's function we deduce the Drude weight which characterizes 
transport even when the description in terms of a mean-free path between two scattering events, 
as in Drude's original theory, breaks down. The dependence of the Drude weight on the boson 
energy uo/tb reveals two opposing tendencies in the correlated transport process: On the one 
hand, transport is assisted by bosons, and therefore D increases with decreasing uio/tb- On 
the other hand, it is the involvement of bosonic excitation in the transport process itself that 
limits its efficiency, and D remains finite in the limit of small wo/tfe- This competition between 
correlations and fluctuations is a general physical phenomenon, which can be further studied in 
the full model [8]. 

The correlated transport process studied here is also important at finite fermion density, when 
the formation of a charge-density wave at half-filling triggers a metal-insulator transition [9j. 
The charge-density wave insulator is peculiar in that it involves few bosons, in contrast to the 
common Peierls scenario, and is caused by strong boson-fermion correlations. A hole doped into 
this insulator propagates through the 6-step-process discussed here, while an additional electron 
can propagate through a more efficient 2-step-process. The simple results obtained here for a 
single fermion are important to understand the more complicated physics in this case. Our 
considerations may also allow for the development of an analytical description of the correlated 
charge-density wave insulator. 
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